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Introduction 


In  this  paper  we  derive  the  probability  distribution  of  customer 
wait  at  a supply  point  following  an  S-1,S  Inventory  policy.  Higa, 
et  al  [1]  derived  the  waiting  time  distribution  for  an  S-1,S  policy  for 
compound  Poisson  demands  and  exponentially  distributed  lead  times. 
Sherbrooke  [3]  investigated  the  distribution  for  the  case  of  constant 
lead  times.  This  paper  derives  similar  results  for  the  pure  Poisson 
case  but  allows  for  an  arbitrary  lead  time  distribution. 

Waiting  Time  Distribution 

Customers  arrive  at  a supply  point  using  an  S-1,S  inventory  policy 
according  to  a Poisson  process.  All  customer  requests  are  for  one  unit, 
and  are  filled  by  a FIFO  priority.  Lead  times  are  lid  with  an  arbitrary 
distribution  F(x). 

We  use  the  analogy  with  the  infinite  server  Poisson  queue  with 
general  ild  service  times.  The  analogy  is  specifically-placing  of  an 
order  : arrival  of  a customer. 

Wa  uaa  Ross'  [2]  pp  17,18  analysis  of  the  infinite  server  Poisson 
queue. 

Definitions: 

x(t)  - number  of  customers  in  service  at  t,  t ^ o 
r(t,x)  ■ number  of  service  completions  in  t to  t + t 

r^(t,x)  - number  of  service  completions  in  t to  t + t of  those  which 
arrived  in  [0,t) 

1 





* v-\-v 


OH 


r2(t,x)  ■ number  of  service  completions  in  t to  t + x of  those  which 


arrived  in  (t,t+x) 

F(z)  - service  time  distribution  function 

Ue  take  t as  the  time  of  a customer's  arrival. 

Now  r(t,x)  - r^t.x)  + r2(t,x)  + «(x) 

where  6(t)  “ 1 if  present  customer's  service  time  is  ^ t 
- 0 otherwise 

The  customer  waits  longer  than  t if  and  only  if 

(1)  x(t)  - r(t,x)  >_  S 

If  the  number  in  service  at  the  customers  arrival  offset  by  returns 
in  t to  t + t is  greater  than  or  equal  to  S , then  there  are  either  back- 
orders which  must  be  filled  before  the  current  order  (x(t)  - r(t,x)  > S) 
or  else  the  next  arriving  item  will  be  used  to  fill  the  current  order 
(x(t)  - r(t,r)  • S).  Since  we  have  Poisson  arrivals,  the  given  condition 
of  an  arrival  at  t has  no  effect  on  the  likelihood  of  past  or  future 
arrivals.  Practically  then,  we  may  deal  with  the  queue  conditions  at  t 
without  regard  to  the  arrival  occurring  there.  Now  (1)  may  be  expressed  as 

x(t)  - r(t,t)  - x(t)  - r^(t,x)  - r2(t,-r)  - S 

But  x(t)  - r^ (t , t ) = Xj^t.x) 

- number  of  arrivals  in  o to  t which  are  still  in  service  at  t + x 


Using  Ross'  method  of  analysis  we  have 


Pr[x1(t,T)  - j|  N(t)  - n]  - (J)  pj  (1-Pj)1 


where 

N(t)  - number  of  arrivals  In  [o,t) 

p^  * probability  that  an  arbitrarily  chosen  arrival  In  o to  t 
Is  in  service  at  t + t 

1 t 1 t+T 

- ± / [1  - F(t+x-z)]dz  - ± / (1-F(z))dz 

O T 

z,  the  time  of  the  arbitrary  arrival  given  N(t)  arrivals,  is  uniform  in  (o,t) 
1 - F(t+x-z)  is  the  probability  that  an  arrival  at  z is  still  in  service  at 
t + x 


Merely  duplicating  Ross'  algebra  we  have  that 

-Xtp.  (Xtpj)-1 

(2)  Fr[x1(t,x)  - j]  - e — j] — 


ie  Poisson  with  parameter  Xtp^  where  X ■ arrival  rate. 
We  also  have  by  similar  reasoning  that 


(3)  r_(t,x)  is  Poisson  with  parameter  Xxp 

A C 


where  p ■ — / F(z)dz  ■ probability 
o 

that  an  arbitrarily  chosen  arrival  in  (t,t+x)  completes 
service  by  t +x. 

Let  Y (t,x)  ■ r2(t,x)  + 6(x) 


* V'-.V 


Then 


Pr [w  > t]  - Pr[x1(t,x)  - Y(t,T)  >_  S] 


Since  x^(t,x)  and  Y(t,x)  are  Independent  we  have 


(4)  Pr [w  > r]  - Z Z Pr[Y(t,x)  - m - k]  Pr[x..  (t,x)  - 
k“S  m-k  1 


On  switching  the  order  of  summation,  and  some  simple  algebra  we  8«t 


Pr [w  > x]  = Z Pr[X]L(t,x)  - S+m]  Z Pr[Y(t,x)  - k] 
m=o  k-o 


Using  the  definition  of  6(x)  we  see  that  Pr[Y(t,x)  - k] 

- Pr[r2(t,x)  - k] [1-F(t) ] + Pr[r2(t,x)  - k-l][F(x)]  for  k > o 


Pr[Y(t,x)  - o]  - Pr[r2(t,x)  - o][l-F(x)] 


Substituting  these  into  (4)  we  have 


Pr[w>x]  - E Prtx^t.x)  - S+m]  E [l-F(x)]Pr[r  (t,x) 
m-o  k-o  Z 


E [F(x)]Pr[r.(tix)  - k - 1] 
k-1  z 


E Pr[Xl(t,x)  - S ' m]  Z [l-F(x)]Pr[r  (t.x)  - k]  + E [F(x>] 
m-o  k-o  k-o 


Pr[r2(t,x)  - k] 

» m-1 

E Pr[x  (t.x)  - S+  m)  [ E Pr[r,(t,x)  - k] ] + [l-F(x)]Pr[r,(t.x)-m] 
m-o  k-o  L 


% V 
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So, 


oo  m-1 

(5)  Pr [»> t ] - Z Pr[x  (t,T)  - S + m]  E Pr[r  (t,x)  - k]  + 
m-1  1 k-o 


Z Pr[x(t,t)  - 8 + m] [1-F(t ) ] Pr[r2(t,x)  - m] 
m-o 


Any  attempt  to  compute  Pr[w  > t]  from  the  above  form  (5)  must  use  an 
approximation.  One  approximation  is  to  use  a finite  sum. 

Let 


N m-1 

(6)  G(t,H)  - E Pr(x..  (t,r)  - S + m]  E Pr[r_(t,T)  - k] 
m-1  A k-o 


N 

+ i Prtx^t.i)  - S + m]  [1-F(T) ]Pr [r2(t,T)  - m] 
m-o 


and 


G(t)  - Pr[w  > T] 


Then  G(t)  - G(T,N)  - E(X,N) 

co  m-1 

- z Pr[x1(t,T)  - s + “I  E Pr[r2(t,t>)  - k] 
tt-N+1  k"° 


+ z Pr(x. (t,x)  - s + m][l-F(T)]Pr[r2(t,T)  - m] 

m-lH-1 


00  oo 

< z Pr[x.Ct,x)  - s + m]  - z Prlx^Ct.-r)  ■ m] 

m-N+1  m-a+N+l 
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Or 


S+N 

E(t,N)  <_  1 - Z Pr [x^(t,t)  **m] 
m**o 

Thus,  we  may  choose  N with  the  assurance  of  achieving  any  error  tolerance. 
Comparison  With  Higa's  Results 

Higa,  et  al  developed  the  customer  waiting  time  distribution  for 
stuttering  Poisson  arrivals  and  exponential  service  lead  times.  Some  re- 
sults for  the  pure  Poisson  case  were  published  in  their  paper.  Following 
is  a comparison  of  their  results  with  this  work  with  exponential  service 
times. 

Arrival  Rate  •>  .01 


Service  Time 

S 

Prob[w  - 
Higa,  et  al 

o] 

Kruse 

Prob [w 
Higa 

< 10] 
Kruse 

5 

1 

.951 

.951 

.999 

.999 

2 

.999 

.999 

1.000 

1.000 

3 

1.000 

1.000 

1.000 

1.000 

10 

1 

.905 

.905 

.986 

.987 

2 

.995 

.995 

1.000 

1.000 

3 

1.000 

1.000 

1.000 

1.000 

20 

1 

.819 

.819 

.928 

.929 

2 

.982 

.982 

.996 

.996 

3 

.999 

.999 

1.000 

1.000 

40 

1 

.670 

.670 

.783 

.785 

2 

.938 

.938 

.968 

.969 

3 

.992 

.992 

.997 

.997 
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The  differences  shown  in  the  above  table,  although  small,  are  real 

differences  between  the  two  models.  They  are  not  due  to  computational 

-4 

approximations,  since  our  resultB  have  been  made  accurate  to  10  . Higa's 

results  are  approximate  because  he  ignores  the  effect  of  what  we  have 
called  Where  we  have  required 

x^t.T)  - r2(t,t)  - 6(t)  >_  S 

for  wait  to  exceed  t,  Higa  requires  only 

- (t)  >.  S 

The  two  models  seem  to  agree  because  of  the  particular  parameter  values 
used  in  the  above  table.  For  these  values,  the  probability  that  ^(t.t) 
is  greater  than  zero  is  always  less  than  .05. 

The  following  table  compares  the  two  models  for  an  arrival  rate  of 
1.0  instead  of  .01.  This  tends  to  make  r2(t,T)  a more  significant  con- 
tributor to  the  probability  distribution  of  wait. 
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ARRIVAL  RATE  -1.0 


Prob  [w-o]  Prob  [w  10] 


Mean  Service  Time 

S 

Higa 

Kruse 

Higa 

Kruse 

5 

1 

.006 

.006 

.806 

.998 

2 

.040 

.040 

.952 

1.000 

3 

.125 

.125 

.991 

1.000 

10 

1 

0.000 

0.000 

.084 

.663 

2 

0.000 

0.000 

.226 

.785 

3 

0.003 

0.003 

.421 

.876 

20 

1 

0.000 

0.000 

0.000 

.005 

5 

0.000 

0.000 

.012 

.085 

10 

0.005 

0.005 

.271 

.501 

15 

.105 

.105 

.790 

.898 

40 

1 

0.000 

0.000 

0.000 

0.000 

25 

0.004 

.004 

.123 

.178 

30 

.043 

.043 

.410 

.491 

40 

.479 

.479 

.933 

.952 

Ignoring  r2(t 

,-c),  then, 

may  in 

some  cases  result 

in  large  errors  in 

the  probability  distribution 

. Note 

that  both  models 

agree  for  G(o) 

since 

' o in  this  case. 
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